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REVISED  AND  ENLARGED  COLLECTION  OF  PLASMA 
PHYSICS  FORMULAS  AND  DATA 


PREFACE 

The  first  edition  of  the  NRL  Plasma  Formulary  was  compiled  and  pub- 
lished in  197^ . By  1975  the  original  printing  of  5000  copies  had  been 
exhausted.  The  time  seemed  ripe  to  prepare  a revised  edition  incorpor- 
ating corrections  of  errors  reported  by  users.  Extensive  new  material 
was  also  added  (some  of  it  already  contained  in  the  companion  manual  to 
the  first  edition,  NRL  Memorandum  Report  Number  269?),  and  a number  of 
changes  made  in  the  interests  of  greater  clarity  and  ease  of  use.  It  is 
hoped  that  the  result  represents  a net  improvement  over  its  predecessor, 
and  that  the  iterations  will  soon  converge  to  a definitive  version. 


Two  forms  of  the  revised  Formulary  were  prepared,  a folding  card 
like  the  197^  edition,  and  a small  booklet.  With  minor  exceptions, 
they  contain  the  same  material.  Their  content  is  reproduced  in  the 
present  memorandum  report,  which  therefore  serves  as  a less  compact 
(but  possibly  more  easily  read)  third  version  of  the  revised  Formulary. 
It  includes  a list  of  errata  for  the  folding  card  and  booklet,  complete 
as  of  the  date  of  publication  of  the  report.  Where  a formula  in  the 


A certain  amount  of  supplementary  material  has  been  included,  particularly 
in  the  Vector  Calculus  section.  The  elaborate  reworking  of  the  Atomic 
and  Molecular  Physics  section  included  in  the  previous  edition  of  this 
report  has  been  abridged,  however.  Interested  readers  seeking  background 
material  on  the  formulas  of  this  section  are  referred  to  the  version 
originally  prepared  by  A.  W.  Ali  for  Memorandum  Report  No.  289S  or  to  a 
forthcoming  report,  "Summary  of  Plasma  Spectrographic  Data,"  by  J.  Davis 
et.  al.,  as  well  as  to  the  numerous  ponderous  monographs  on  the  subject. 

Once  again,  many  of  my  colleagues  in  NRL's  Plasma  Physics  Division 
helped  with  the  Formulary'  by  suggesting  changes  and  corrections  and  by 
proof  reading.  In  addition  to  those  persons  previously  acknowledged, 
particular  thanks  for  specific  contributions  are  due  to  Dr.  Verne  Jacobs 
of  Science  Applications,  Inc.  (Atomic  Physics  and  Radiation'' , and  Drs . 
Niels  Winsor  and  Stephen  Bodner  of  NRL  (Collisions  and  Transport  and 
Lasers,  respectively’).  My  special  thanks  go  to  Ms.  Jean  Monnette  for 
managing  to  compose  and  lay  out  the  Formulary  while  I was  making  con- 
tinual changes  in  it. 


ERRATA 

March  1,  1977 


(Foldout) 

Panel 

2 


(Booklet) 

Page 


Erratum  and  Correction 
Eq . (11):  replace  _Bx7xA  by  _Bx(7xA' 


Eq.  ( 17-19) > (27):  See  p.  6 


5 


lk 


16 


19 


7 Under  "Components  of  (A  • 7)  B": 

oEz  cB 

replace  by  

c oz 

z 

10  Under  "Rationalized  mks" 
replace  "mho"  by  "siemens" 

11  Under  "Rationalized  mks" 
replace  "weber/m2"  by  "tesla" 

12  Under  "Rationalized  mks" 
replace  "newton/m2"  by "pascal" 

jb  Ion  gyrofrequency , F region: 

replace  220  by  JOG 

21  Def.  of  L:  "horixontal"  should  be  "horizontal" 

33  In  Weakly  Ionized  Plasmas,  Eq . for  collision 

frequency  v : T.  should  be  multiplied  by  k. 

Ct  oc 

Eq . ( c':  interchange  11.°  and  2 .u ; last  Eq . : 

replace  A E by  A / E 

bh  Eq . ( 1:  ':  should  read  = . 71'  x 1 1 ‘ T 

e 

cm  /sec 


p 


NUMERICAL  AND  ALGEBRAIC 

Euler-Mascheroni  constant^  v = 0.57722 
Gamma  Function  T (x  + 1)  = x F (x) 

r(i)  = 5.5665 
r(|)  = 4.5906 

r(  = ) = 5-6256 

4 

r(|)  = 2.57S9 
r(|)  = 2.21S2 
r(i)  = 1.7725  = tt- 
r(|)  = 1.4892 
r(|)  = 1.5541 
r(f)  = 1.2254 
r(|)  = i.i542 
r(-J  = 1.1288 
r(i)  = 1.0 

Binomial  Theorem  (good  for  |x|  < I or  Ci  = positive  integer': 


(1  + x)a  = ! (?)xk  , 1 + a x ♦ ^ (P~g)  x-  + 

k=0  2 ! 5 ! 


Here 


( 


= r(  r+l'i  = q(q-  1' . . ■ (q-k+1' 
k . ' k!F( '+1- k'  k! 


is  the  binomial  coefficient. 

A verv  general  relation  involving  the  binomial  coefficients  is  the 

[2] 

Rothe-Hagen  identity  (good  for  all  complex  x,v,z  except  when  singular': 


r x / x+kz  \ y /y+(n-k)z\ 

~ x+kz  \ k / v+(  n-  k ' z \ n-k  / 

n -1  v 


x4y 

x+-hiz 


Manv  formulas  involving  binomial  coefficients  are  special  cases  of  the 
Rothe-Hagen  identity.  For  example,  with  z = 0, 


JoOO(-)  - 


with  x = y = n,  z = 1, 


£ ( n+k-  1 1 ! (2n-k-l'>  I 

k=  0 k!  (n-kl! 


2 n 1 ( pn)  ! 
~ n^  (2n) ! 


with  x = - y = — , z = — , n odd, 

o r, 


k=0  n+2k  n-2k 


/ n + 2 k\  /n  - 2 k\ 


ki  \ 

/ ' 


i 


VECTOR  CALCULUS 


Notation:  f,  g,  etc.,  are  scalars;  A,  B,  etc.  are  vectors;  T is  a tensor 

(1)  A-B  x C = A x B-C  = 3-C  x A = 3 x C-A  = C-A  x B = C x A-B 

(2)  A x (B  x C)  = (A-C)  3 - (A-B)  C 

(3)  A x (B  x C)  + B x (c  x A)  + C x (A  x B)  = 0 

(t)  (A  x B)  • (C  x D)  = (A-C)  (B-D)  - (A-D)  (B-C) 

(5)  (A  x B)  x (C  x D)  = (A  x B-D)  C - (A  x B-C)  D 

(6)  v(fg)  =7  (gf)  sf  78+g7f 

(7)  5-(fA)  = f 2-A  + A.r  f 

(S)  7 x (fA)  = f 7 x A + 7f  x A 
(9)  7. (A  x B)  = B^.7  x A - A. 7 x B 

(10)  7 x (A  x 3)  = A (7.3)  - 3 (7 .A)  + (3.7>A  - (A.V)  B 

'(11)  7 (a*3)  = A x (7  x B)  + B x( 7 x A)+  (B-7)  A + (A-7)  B 

(12)  72  f = 7-7  £ 

(13)  72  A = v(7-A)  - 7 x 7 x A 
(It)  7 x £f  = 0 

(15)  7-7  x A = 0 

If  e , e , eo  are  orthonormal  unit  vectors,  a second-order  tensor  T 

“1  “2  "3  - 

can  be  written  in  the  dyadic  form 

(lo)  T = Z T. .e.e. 

In  cartesian  coordinates  the  divergence  of  a tensor  is  a vector  with 

components  given  by 

(17)  (7  • T).  = Z (3T  . ./Sx.) 

- : 1 j J3-  J 

In  general 


(It'1  7-(AB)  a (7  • A)B  + (A  * 7)B 

(19)  7-(fT)  a 7f  • T + f 7 . T 


Instead  of  (17),  an  alternative  definition  '■  was  employed  in  the 
pocket  versions  of  the  Formulary: 

(17’''  (7  • T)  = I (oT. ,/ox.) 

- = i j ij  J 

Since  we  vnrite  the  operator  7 to  the  left,  this  choice,  involving 
contraction  with  respect  to  the  second  index  instead  of  the  first, 
appears  less  "natural."  When  it  is  adopted,  ( 18)  and  (19) 
become 

( I- 1 ) 7 • (AB)  = A 7 • B + B • 7 A 
(19')  7 • ( f j)  = T • 7f  + f7  • T 

But  Gauss'  theorem  for  tensors,  formula  (27"',  becomes 

(27')  ? 7 • T dV  = / T • dS  , 

V S 

which  to  the  unsuspicious  eye  is  equally  "natural"  looking.  For 
symmetric  tensors,  the  primed  and  unprimed  formulas  are  identical. 

In  plasma  physics,  fortunately,  most  of  the  tensors  we  deal  with  are 
symmetric.  Nevertheless,  the  choice  fl?1'1  is  a potential  source  of 
confusion  and  its  use  in  the  other  editions  is  regrettable.  Here  I 
have  abandoned  it.  For  uniformity,  all  of  the  integral  formulas 
(25-3*0  have  been  written  with  differentials  at  the  left,  although 
this  results  in  a significant  change  only  in  (27).  The  tensor 
divergence  7 . T in  cylindrical  and  spherical  coordinates  (p.  12  and 
p.  It,  respective  ly)  now  appears  with  indices  transposed. 


— — - -is mrmmm. 


Let  r = ix  + jy  + kz  be  the  radius  vector  of  magnitude  r,  from  the 
origin  to  the  point  x,y,z.  Then 

( 20 ) V«r  = 3 

(21)  7 x r = C 

(22)  7 r = r/r 

(25)  7 (1/r)  = -r/r3 

(24)  7.  (r/r3)  = 4ri(r) 

If  7 is  a volume  enclosed  by  a surface  S and  dS  = n dS  where  n is  the 
unit  normal  outward  from  V, 

(25)  J*  dV  7 f = J* dS  f 


(26)  J dV  7 • A = J*  dS  • A 


(27)  P dV  7 *T  = f dS  • T 


(28)  f dV  7 x A = P dS  x A 

O — — J — — 


(29)  J dV  (f  fg-  g]f  f)  = JdS  • ( f 7 g - g7fl 


(30)  Jdv(A.7x7xB-f,x7xA)  = P dS-(Bx7xA-Ax7x 

XT  S 


If  S is  an  open  surface  bounded  by  the  contour  C of  which  the  line 
element  is  dl  , 


P3  l 


Cylindrical  Coordinates 


r 


L 


Spherical  Coordinates 


Divergence 


1 3 1 

7 -A — (r2A  ) + 


2 -n  r x a 

r 3r  r sin  9 ?9  ’ 


3 1 

— ( A sin  9)  + 


cA 


r sin  6 5^ 


Gradient 


( Vf ) = — ; ( 7f ) = I M 

- r „ 1 _ ' 9 


Sr 


r 59 


= 


1 Sf 
r sin  9 5«$ 


Curl 


1 3 1 aa 

(7  x A)_  = — (A,  sin  9) — 1 


r ' a 

r sin  9 59  p 


r sin  9 5^ 


, 1 3A  1 5 

(7  x A)  = — - ~ — (rA,) 


9 , '•‘‘"■A1 

r sm  9 54  r Sr  0 


1 3 


, - 1 5A 

(7  x (rA  ) £ 

r 3r  r 59 


Laplac  ian 


T^f  = — — 


5 f 


i r — + 


3r  \ 5r 


Af  \ 

sin  9 — ! + 


1 3 f 


r sin  9 SS  \ 59  / r2  sin2  9 54s 
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Laplacian  of  a vector 


(72A'I  = 72A 


2A  2 3AQ  2A  cot  9 

r y tr 


BA 


r r2  r2  B9 


r~  sin  9 B«$ 

2 cos  9 BA, 


2 BA  A.  ^ '-'-'o  « 

(^A'l  = ^A.  + — — — 

- ® ® r2  B9  r2  sin^  9 sin2  9 ai 


Components  of  (A-V'iB 


- - - r r 


3r 


o 

d 

BA 

r 

sin’  9 

T 

r~  sin  9 

u 

BB 

r 

t ^ 

3B 

r 

r 

B9 

r sin  9 

3^ 

A 

aB9 

!!i 

r 

o© 

+ 

r sin  9 

3<4 

A 

+ ^ 

r 

B9 

r sin  6 

B«S 

iA  2 cos  9 BA. 


-2  „ ,•  _2 


9 9 a6  6 


It 


General  Curvilinear  Coordinates L 1 J 

In  any  orthogonal  curvilinear  coordinate  system,  all  of  the  opera- 
tions involving  7 can  be  defined  in  terms  of  the  scale  factors  Ik. 

These  are  related  to  the  differential  arc  length  dX  by 


(35) 


dX^  = dx  + dv“  + dzc 


h2d?£  + h2d?2  + h2d? 

1 "1  2 '2  3 


where  > designate  the  curvilinear  coordinates.  For  example, 

t '1  ”2  *3  ) 

the  scale  factors  associated  with  the  spherical  coordinates  j r,  0,  z j- 
are  1,  r and  r sin  S,  respectively. 

In  orthogonal  curvilinear  coordinate  systems  a vector  A can  be 
written  as 


(5o) 


A = eA  +eA  +eA  , 

'll  "22  -3  3 


where  the  | e^  | form  a right-handed  set  of  unit  vectors  normal  to  the 
surfaces  of  constant  Then  the  divergence  of  A is 


(37) 


7 -A 


where  H = h h h . The  gradient  of  a scalar  f is  given  by 

12  3 


(36) 


1 


From  Eqs  • ( 1— ^ snd  (y?)9  il  follows  that  the  Laplacian  of  f 

(39^  72f  = - V JL  (l L 1L  \ 

H Y ae.  I h!  a?.  ; 

i - i \ i =i  ' 

The  curl  of  A is 


(V  X A).  = y e.  ..  1 3 fu  ^ \ 

- - i ijk (h  A ) , 

j,k  h.h,  3?.  K K 

J k -j 


are  not 


where  equals  + 1 if  |ijk|  is  even  ^i.e.,  an  even  permutation  of 

{l23j),  equals  - 1 if  |ijkl  is  odd,  and  vanishes  if  i,j  and  k are  no I 
all  distinct. 

To  obtain  expressions  for  the  remaining  operations,  it  is  useful 

to  introduce  the  Christoffel  svmbols  F^,  , defined  bv 

jk 

, 1 ! 3h . . 3h,  , : 

<*»  r k — ■ 6 - -S  6.1  ! , 

*J  k _ a?k  J 

where  5.  equals  1 if  i = j and  vanishes  otherwise.  Fi  vanishes  if  i 
J jk  J 

and  k both  differ  from  i or  if  j = k.  In  addition. 


(42) 

r*  + r* . = 0 ; 

jk  kj 

(43) 

r1,  + rl . + rk . = 0 ; 

jk  ki  ij 

and 

(441 

Y'  r-i  _ 1 3 H ' 

^ ‘ Ik  ~ „ , 

1 H r:k  hk 
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and  less  in 


There  are  at  most  six  independent  non- zero  componen 
most  coordinate  systems  of  physical  interest.  The 


ts  r ..  , 
jk’ 


curvature 


tensor , 


(1+5) 


Ri  jk£ 


i ar; 


k i 


h.  3c. 
J J 


+ 


rn  rn  _ r-i  r*-  + 1 

I i j ' k£  ‘ nk  n.C  ‘ kn  Hr.  j 
n 


always  vanishes  in  a flat  space,  i.e.,  one  which  can  be  described  by 
Cartesian  coordinates.  (This  is  the  only  kind  of  space  we  usually  deal 
with  in  plasma  physics.) 

In  terms  of  the  F's,  the  gradient  of  a vector  A is  a tensor  T 
given  by 


(43) 


T.  . 
Ji 


(V  A) 


J + 


3A . 


h.3£. 
J 'J 


Contracting  this  expression  (i.e.,  multiplying  by  5^  and  summing  over 
i and  j)  recovers  Eq . (57)-  The  directional  derivative  is 


7)B 


A.  3B.  A.B.  3h. 
_J.  _i  + -LJ.  i 

h.  3c.  h.h.  3r . 
J 'J  1 J 'J 


A .B  . 3h  . 


The  gradient  of  a (second-order''  tensor  is  a third-order  tensor  given 

by 


l 


9T  . 


(43) 


Cl  V 


V!k 


Z 


rk. 

a 31 


+ rk  T \ 

+ ' ULZi  | 


The  divergence  of  the  tensor  T is  obtained  by  contracting  Eq . - with 

respect  to  the  first  two  indices: 


If  T is  defined  according  to  Eq.  , 7 . T is  the  Laplacian  of  the 

vector  A: 


<^'1 Z z 

j j 


1 9 

H 95. 
'J 


r 

I H 

' h? 

L J 


oA . 

1 

a?. 

2 j 


(50> 


h.A. 

L 1 


95  . 

J 


9h. 


hThT  9?. 
t ] J j 


9h.  9 

— - (h  .A  .) 

2 2 


hfh “f  9?.  95. 
t J J 1 


9h . 9 h .A . 9 19 

— -1  (h.A.")  + — ^ ^ ~ 

h.h3  9?.  9?.  J J h 9?  H 9? 

i j 1 j 1 1 _ j 


Evidently  the  first  term  on  the  RHS  is  v A, 


to  A . 
1 


Eq . (39''  • Fortunately,  Eq . ("O')  s 
most  coordinate  systems  of  interest. 


, the  scalar  Laplacian 
implifies  considerably 


applied 

for 
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In  the  following  section,  the  results  of  Eqs.  (55)  - (50)  are  applied 


to  a coordinate  svstem  occasionally  encountered  in  plasma  physics 
problems . 

Toroidal  Coordinates 

An  arbitrary  point  in  5-space  can  be  uniquely  identified  by  the 
coordinates  ■*  o , A,  'Jr  j.  , where  o is  the  distance  between  the  given  point 
and  a fixed  toroidal  axis  of  radius  R (the  major  radius),  A is  the 
azimuthal  or  toroidal  angle  (i.e.,  Che  one  sweeping  around  the  major 
circumference'1,  and  jr  is  the  poloidal  angle.  (Note  that  this  system 
differs  from  the  "toroidal"  coordinates  in  Ref.  5,  p.  ■:-:6f  and  from 
the  flux  surface  coordinates  commonly  used  in  connection  with  toroidal 
devices.''  The  values  assumed  by  these  coordinates  are  restricted 
according  toOio<cc,  0 £ A < 2t,  and  0 ^ jr  < 2tt.  If  j r , A,  z | are 
the  cylindrical  coordinates  of  the  point  in  question,  where  z is  normal 
to  the  plane  of  the  torus,  then  A is  the  same  in  both  systems  and 

R + o cos  t = r, 
o sin  i = z . 

Hence  the  square  of  the  differential  are  length  is 
dir  = dr  + r^dt^"  + dz~ 

- do  + (R  + o cos  5)2d«4“  + o‘~ dy“  . 

Consequently  the  scale  factors  are  1,  R + o cos  i,and  c,  respectively. 


From  Eq . (37'),  the  divergence  of  A is 

7 'A  = j — — - o(  R + o cos  •;')A  + 

R + 0 cos  i><  i n Ac 


i h 


3A,  1 o r , 

H — — (R  + o cos  %)A'i  ] 

n A if:  * J ! 


3<5  DO 


From  Eq . (35'',  the  gradient  of  f has  components 

(7f)  ; 


R + o cos  »!r  06 


(7fV 


From  Eq . (39"),  the  Laplacian  of  f is 


7£f  = j — — o(R  + o cos  it)  — 

R + d cos  'll  I o co  3o 


R + o cos  3<5'  O’1  S'l 


M J I 


‘ 1 


From  Eq.  (^-0),  the  curl  of  A is  given  by 


1 (BA,  IB 

(7  x A).  = — - - - — 

( B 6 o 3 


- 6 


(7  X A) , = 


R -t 

- 0 COS  'i 

1 1 

l BA  B 

0 

o 1 

1 B*!r  Bi 

i 

(R  + 3 cos  ill' A, 


J| 


•’I 


' R + o cos 


•Jf  I 


i 9 

1 

> 
1 

) 

(R  + o cos  ir^A/  - 

L 6 u- 

I Bo 

The  only  nonvanishing  Christoffel  symbols  are 


= _ p2  = COS 

* 2i  i; 


R + o cos  •> 


r2  = - r2 

23  32 


- sin  •!; 

R + p cos  t 


= _ p3 

*31  *13 


BB_  Aj5  SB. 


Then  from  Eq.  (^-7^, 

[(A  - ZBb1  =ao^+ ^ 

L Jo  Be  R + o cos  'if 

A 3B  A ,B  , cos  v A . B . 
o o 0 v 


R + o cos  v 


L ‘ 

- -u 

p 

3o  R 

+ 0 cos 

u 

\ 

A ,B 
, 0 c 

cos  'i  A^B,. 

sin  t 

o 

3-i- 

R + 3 

cos  'i  R + 

C COS  *; 

A . SB  , A ,B  A ,B  , sin  * 

■It  £ + 'it  o + p g 

o 31-  o R + o cos  >!f 

From  Eq . (^9^  > 

(<?  ■ t)  = ’ — — o(R  + o cos  -iO T 

R + c cos  '1  I o 3c 


b R + c cos  *!;  ! 


Finally, 

( 


-T^  1 
+ — — + - 


o si- 


(R  + o cos  t)T 


I 


+ cos  - sin  ; 


* 


(V  . T'1  = 


I 1 


R + 3 cos 


i I 3 


o(R  + 3 


oT,  1 3 

+ — ^ + 


bi  a c'h 


"R  + o cos  ’lOT 


+ sin  *T , , . + — T 

66  1 to 

’ o 


from  Eq . ( 5 CO  , 


=A'  = -^A  - A 


-1  + 


cos' 


(R  + 3 cos  t'*2  J 


COS  7 


^ 2 oA, 


PHYSICAL  CONSTANTS  > I 

CSS 


Physical  Quantity 

Symbol 

Va  1 ue 

Units 

Boltzmann  constant 

1 

k 

1.5807  X icr1: 

erg/deg  ( I<) 

Elementary  charge 

e 

4.8032  x lcr10 

statcoulomb 

Electron  mass 

m 

9 . 1C95  x lc-^ r 

g 

e 

Proton  mass 

mP 

1.6726  x lc-24 

g 

Gravitational  constant  , 

G 

6.6720  x icr5 

dyne-crrr/g“ 

Planck  constant 

h 

6.6262  x icr2T 

erg- sec 

h = h/2rr 

1.0546  x lcr-7 

erg- sec 

Speed  of  light  in  vacuum 

c 

2.9979  X 1010 

cm/ sec 

METRIC  PREFIXES ^ 


Multiple 

Prefix 

Symbol 

Mul tip le 

Prefix 

Symbol 

10'1 

deci 

d 

10 

deca 

da 

-o 

10 

centi 

c 

102 

hecto 

h 

10-3 

milli 

m 

103 

kilo 

k 

10"  ' 

micro 

U 

10J 

mega 

M 

10’9 

nano 

n 

109 

giga 

G 

1C  12 

pico 

p 

c\j 

H 

o 

r— < 

tera 

T 

10'15 

femto 

f 

1015 

peta 

P 

10-1  : 

atto 

a 

10ls 

exa 

E 

ELECTRICITY  AND  MAGNETISM 


In  the  following,  c = conductivity,  f - j-A  - radiation  frequency, 

'f  = l:/u  and  '<  = e/s  . Where  subscripts  are  used,  1 denotes  a con- 
ni  o e o 

ducting  medium  and  2 a propagating  lossless  dielectric  medium.  All  units 
are  mks  unless  otherwise  specified. 


Permittivity  of  free  space 


s = 2.6542  x 1CT12  farad/m. 
o 


Permeab i 1 i tv  of  free  space 


10" 7 = 1.2566  x 10"'-  henrv/m. 


Resistance  of  free  space 


R = (U  /e  7‘/s 
o 00 


= 57-. 75  onm 


Relaxation  time  in  a lossy 
medium 


s /c  sec 


1/2  _i  /2 

Skin  depth  in  a lossy  medium  s = (2/titicr^  = (-fur)  m 


Wave  impedance  in  lossy  medium 


= [ u/(  e -hie/ id)  J 1^2 


Reflection  coefficient  at 
conducting  surface  (good 
onlv  for  R s 1) 1 ^ 


R = 1 - 4.22  x 10"*(fC  '<  fa) 

mi  er 


Field  at  distance  r from 


straight  wire  carrying 


current  I 


B„  = ul/5— r tesla 

= 5.2  i/r  gauss  (r  in  cm) 


Electromagnetic  energy  in 


vo lume  V 


If/ 

W = - dv( H • B + E-D7 


Povnting's  theorem 


— 4-yVdS  = --jt dV J • E 


where  S is  the  closed  surface  bounding  V and  the  Poynting  vector  energy 


flux  across  S is  N 3 E X H 


MAXWELL'S  EQUATIONS 


Name 

Rationalized  mks 

Gaussian 

Faraday's  law 

BB 

1 '-B 

7 x E = - — 

7 

x E = - 

0 t 

' 

c 

Ampere ' s law 

?D 

7 x H = — 4 J 

7 

1 5D  4- 

x H = — — n.  + — — j 

' * at  ' 

c at  c - 

Poisson's  equation 

7-D  = o 

7-D  = 

£ Absence  of 

7 • B = 0 

O 

II 

CO  I 

t>  1 

magnetic  monopole  ] 

Lorentz  force  on 

q £ E + v x B 1 

q 

[ ? +i Y X 3 ] 

charge  q 

Constitutive  relations 

D = e E 

D = e E 

B = UH 

B = uH 

In  a plasma,  u *5  ^ (Gaussian  units:  u » 1).  The  permittivity  satisfies 
s (Gaussian:  € ss  1}  provided  all  charge  is  regarded  as  free.  Using 

the  drift  approximation  v = E x B/B""  to  calculate  polarization  charge 
density  given  rise  to  a dielectric  constant  K = e/e  = 1 + J>6  r x lcF  o/Bc 
(mks)  = 1 +4— cc2/B‘  (Gaussian),  where  o is  the  mass  density. 


37 


ELECTROMAGNETIC  FREQUENCY /WAVELENGTH  SALTS 


^The  boundary  between  ULF  and  ELF  is  variously  defined. 

^The  SHF  (microwave)  band  is  further  subdivided  approximately 

rsl 

as  follows r1- 


Sand 

Frequency  Range  (GHz) 

Wavelength 

Range  ( cm) 

Designation 

Lower 

Upper 

Lower 

Upper 

S 

2.5 

3-95 

7.5 

11.5 

G 

3-95 

5.85 

5.1 

7.6 

J 

5.3 

8.2 

3-7 

5-7 

H 

7.05 

10.0 

3.0 

4=“ 

f 0 

vn 

X 

8.2 

12.4 

2 . “p 

3-7 

M 

10.0 

15.0 

2.0 

3-0 

P 

12.4 

18.0 

1.67 

2.4 

K 

18.0 

26.5 

1.1 

1.67 

R 

26.5 

40.0 

.75 

1.1 

AC  CIRCUITS 


(c) 


L'nderdamped,  £ < 0 


qt  = 


- - , "V 

hi  -(Y-,q  + I ) sin  ai  c r q cos  s t | e 

R o o i ~ o 


wnere 


ci)  = U)  ( 1 - R2cAL  )Z/2  , Cl  = (LC)-1/2 
1 o \ ' o 


The  quality  of  the  circuit  is  Q * cju^L/R ; cu  is  the  resonant  frequency. 
At  u)  * cu  , 2 = R.  Instability  results  when  L,  R,  C are  not  all  of  the 
same  sign. 


DIMENSIONLESS  NUMBERS  OF  FLUID  MECHANICS 


Hundreds  of  dimensionless  groups  have  been  defined,  particularly  in 

f - o i 

chemical  engineering  and  fluid  dynamics. ” ' The  present  collection 

includes  only  a small  subset  of  the  latter  category. 

Notation  is  as  follows: 


c 

c 


g 

H 

k = v.cc 

P 

L 

n 

r 

R 

a 


speed  of  sound 

specific  heat  at  constant  pressure 
(units  cnr/sec~  - deg) 

gravitational  acceleration 

vertical  length  scale 

thermal  conductivity  (units  g-cm/sec  '-deg) 

•SI* 

horizontal  length  scale 
frequency  of  solid-body  oscillation 
radius  of  curvature  of  a pipe  or  channel 
radius  of  pipe 

Newton's  law  coefficient  of  heat  transfer 
(units  g/sec  -deg) 

volumetric  expansion  coefficient,  dV/V  = .^dT 

applied  temperature  difference 

thermal  diffusivitv  (units  cm/sec) 

collisional  mean  free  path 

bulk  viscosity 

kinematic  viscosity 

rotational  velocity 


mass  density 


Same 


Symbol 


Definition 


S igni finance 


Burgers 

B 

(3gAT)*/2nH* 



effect  of  stratification/effect  of 
rotation 

Dean 

D 

i- 

(UR/v)(R/r> 

transverse  circulation  due  to 
curvature/ longitudinal  flow  speed 

Eckert 

E 

irVc  AT 

P 

adiabatic  temperature  increase/applied 
temperature  difference 

Ekman 

Ek 

v/:m2 

viscous  forces/coriolis  forces 

Froude 

Fr 

^L/g 

centrifugal  forces/gravitat ional 
forces 

Grashof 

G 

gpMTL 7 v2 

effects  of  bouvar.cy/viscous  effects 

Knudsen 

Kn 

X/L 

hydrodynamic  tirre  scale/collisional 
time  scale 

Mach 

M 

U/c 

directed  motion/thermal  motion 

Nusse It 

N 

-OfL/k 

total  thermal  transfer/ thermal 
conduction 

Pec  Let 

Pe 

UL/h 

thermal  diffusion  time  scale/free 
streaming  time  scale 

Prandtl 

Pr 

v/n 

momentum  diffusion  rate/heat 
diffusion  rate 

Rayleigh 

Ra 

rgiTH  7vk 

bouyancy  forces/diffusive  forces 

Reynolds 

Re 

GL/v 

inertial  forces/viscous  forces 

Richardson 

Ri 

SgATH/U2 

effects  of  bouvancv/ef fects  of 
vertical  shear 

Rossby 

Ro 

u/r.L 

inertial  forces/coriolis  forces 

Stanton 

s 

X 

o/Pc  U 

p 

thermal  convection/compress ional  heating  j 

Strouhal 

s 

nL/U 

oscillatorv  mot  ion  Mi rected  motion  ( 

r 


SHOCKS 


Ideal  gas  jump  conditions  for  density  o,  velocity  u,  pressure  p, 
temperature  T: 


(y  + 11  M2 


p 2 vM2  - v + 1 


u ( y - ll  M2  + 2 p 

2 11 


V + 1 


T 2 vM2-  v + 1 ( v - llM2  + 


(y  + D2m2 


where  1,  2 denote  quantities  in  front  of  and  behind  the  shock  respectivel\ 

in  the  shock  frame;  M = u /c  = Mach  number  in  region  1,  c 2 = vkT/o; 

1 1 si  s' 

V = ratio  of  specific  heats  = ' f + 2)/f,  where  f is  the  number  of  degrees 
of  freedom.  Entropy  change  across  the  shock  is 


£S  = S - S * c £r.  — / 

2 1 v _ 


pP 

P \ 0 

L t ' o 


where  c^  is  the  specific  heat  at  constant  volume.  For  weak  shocks, 


c 11  (M2  - 1): 

v 5 (v  + l)2  1 


Evidently  in  the  limit  M -•  p /o  -*  (y  + l1/(' 


while  the  temperature  and  pressure  jumps  become  infinite. 


PLASMA.  DISPERSION  FUNCTION 


fl  1J 


Definition 


Z(0  = J~ dt  exp( - t2V(t  - C ) (lm£  > 0) 
- 00 

iC 

= 2 i exp(  - C !2) dt  exp(-  t2")  . 


Physically  ' = x + iv  is  the  ratio  of  phase  to  thermal  velocitv. 


Differential  equation 


dZ 
d C 


= -2[l  + CZ],  Z(C^  = iTT*  ; 


£2  + 2 r 12  + 

d;2  dr 


2 Z = 0 . 


Real  argument  ( v = O' 


Z(x'  = i77"'  exp(-x" ' - 2 x Y(x)  , 


y(x)  = x_1  exp( -x2'' J*dt  exp(t“). 


Imaginary  argument 


Z(  iv 


ir  exp y'  [l  - erf(v)’j  = i~”  exp(y  'erfc'v' 


Ls 


Power  series  small  argument 


Z ( 0 = irr~  exp (-  C2  ' - 2 £ | 1 - 2 + 4 C*/15 

L 


k . 


FUNDAMENTAL  PLASMA  PARAMETERS 


All  quantities  are  in  Gaussian  units  except  temperature  ' Tg , T^,  T' 


expressed  in  eV  and  ion  mass  (nu)  expressed  in  units  of  proton  mass. 


U = m^/m^;  2 is  charge  state;  k is  Boltzmann's  constant;  K is  wavenumber. 


Frequencies 

electron  gvrofrequency 


ion  gvrofrequency 


hybrid  gvrofrequency 


electron  plasma  frequency 


ion  plasma  frequency 


f = (u  /2r  = 2 .tO  x 10;  B Hz 
ce  ce 


uu  = eB/m  c = 1 .Jo  x 10'  B rad/sec 
ce  e 


f . = 0)  . /2rr  = 1.52  x 10  Zu-1  B Hz 
ci  ci 


0)  . = eB/m.c  = 9.58  x 10"  Zj”“  B rad/sec 
ci  1 


fm  = tlh/2_  = -53  x 104  z*  u“* 


U ? B Hz 


i'TU  = \U)  x-  . 

LH  ce  ci 


, 10  x 10~  Z ' u * B rad/ s 


f = uu  /£"*  = 5.98  x ICt  n * Hz 
pe  pe  e 


pe 


( Vrn  e /m 


5 . x 10"*  n r rad/ sec 
e 


f . = uu  ,/2~  = 2.10  x 10=  Zu“*  n ~ Hz 


pi  pi 


0)  . = (4rrn.  Z‘  e /’m. 
pi  11 


-1  _ i i 

l.;2  x 10'  Zu  n.  rad/sec 

l 


electron  trapping  rate 


vT  = / ^ = 7-26  x 10e  sec 

Te  l «.  / 


ion  trapping  rate 


Velocities 


(?)*- 


, -r  i-  i 

b 9 x 1C7  K E ' ~ sec 


electron  thermal  velocity 


ion  thermal  velocity 


ion  sound  velocity 


Alfven  velocity 


\ 1 

v_  = (kT  /m  ) " = 4.19  x 107  T c cm/sec 
Te  e e e 

— i t 

v ^ = (kT^/nOs  = 9 - 79  x 105  U ■ T.'  cm/sec 


c = (kT  /m.)  ' = 9*79  x 10s  u - T ^ cm/sec 
s e i e 


v^  = B / ( 4-TTn  .tiu  ) « 

_ l _ l 

= 2.18  x 1011  Li  r n.  ? B cm/sec 

l 


With  this  definition  of  thermal  velocity,  the  Debye  length  and  plasma 
frequency  for  a species  j are  connected  by 


V»J  ■ VT)  ■ 


where  is  obtained  by  substituting  T.  for  T in  the  definition  below. 
Dj  7 ° J 

Some  individuals,  however,  prefer  the  definition 

v.^  = (2  kT7nu)^  , 

according  to  which  the  electron  and  ion  thermal  velocities  and  ion  sound 
velocity  become 


v'  = (2  kT  /m  )E  = 5-93  x 107  T s cm/sec; 
Te  e e e 


v . = (2  kT./m.)?  = 1.3-  x 10"  (T./u)?  ctn/sec; 
Ti  1. 


c = (2  kT  /m.y 
s e 1' 


= 1.3c  x 106  (T  /ur  cm/s 


The  Debye  length  and  gyroradii  may  or  may  not  be  redefined  accordingly. 


Lengths 


electron  deBroglie  length 


a = */(m  kT  )"  3 2.7o  x 10"s  T "5  cm 
e e e 


classical  distance-  of 
minimum  approach 


a2/kr=  1 M x icr7  r1 


electron  gyroradius 


r 3 v„  /uu  3 2.36  T c B-1  cm 
e Te  ce  e 


ion  gyroradius 


plasma  skin  depth 


Debye  length 


r.  = v /ou  . 

1 Ti  ci 

1 1 1 

3 1.02  x 102  U*  z'7  T.s  B_1  cm 

1 

cr  _ 

c/a)  3 5.31  x 103  n ' cm 

pe 

1 1.  I 

X 3 (kr/l+TTne2)  ' 3 7. 1+3  x 1 02  T n'? 


Dimens  ionless 


electron/proton  mass  ratio")  ' 


(m  /m  ) ' 3 2.33  x 10'"  = 1/42 .9 

e p 


number  of  particles  in  Debye 
sphere 


— n\®  3 1.72  x 1CP  TS/2  n-? 

7 ^ 


Alfven  velocity/speed  of  light 

v / c = 7.28  if*  n.“*  B 
A * 1 

magnetic/ion  rest  energy  ratio 

B2/ern.m.c  = 26.5  lT1  n."1  B2 
11  1 

electron  plasma/gyrofrequency 

UL>  /x  = 5.21  x 1CT3  n f B-1 

pe  ce  e 

ratio 

ion  plasma/gyrofrequency  ratio 

x . /x  . = 0.157  if  n.~  B-1 

pi  Cl  1 

thermal /magnetic  energy  ratio 

0 a SrrnT/B2  = ^.03  x 1CT11  nT  B-2 

Miscellaneous 

Bohm  diffusion  coefficient 

ckT  „ „ 

D - = o.25  X 10  TB  cnT  sec 

B 16  eB 

Spitzer  resistivity 

T = 1.15  x 10"14  Z in  A T" 3 7/2  sec 

( transverse) 

= 1.03  X 10"2  Z in A T"3'2  ohm- cm 

The  anomalous  collision  rate  due  to  low  frequency  ion  sound  turbulence 


v*  * x W/kT  =■  5.54  x 104  n_1/2  W/kT  sec"1, 

pe 

where  W is  the  total  energy  of  waves  with  x/ K < VT^* 

Magnetic  pressure  is  given  by 

P = 5-96  V 10  - B‘  dyne/cm2, 

when  P anu  B are  expressed  in  cgs  units.  Thus  the  pressure  equivalent 
of  10  kC  is 

Pf 10  kC)  = 5 .96  x 10  dyne/cm  = T .°5  atm. 
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APPROXIMATE  MAGNITUDES  IN  SOME  TYPICAL  PLASMAS 


'Oflio  n Icin'1) 


Fig.  1 — Information  of  a similar  nature  is  depicted  graphically"* 


COLLISIONS  AND  TRANSPORT 


Temperatures  are  in  eV ; the  corresponding  value  of  Boltzmann's 
constant  is  k = l.oO  X 10  ~~~  erg/eV;  masses  |J,  a'  are  in  units  of  the 
proton  mass;  e,  = Ze  is  the  charge  of  species  Of.  All  other  units  are 
cgs  except  where  noted.  ^ 

Relaxation  rates A 

Rates  are  associated  with  four  relaxation  processes  arising  from  the 
interaction  of  test  particles  (labeled  0!)  streaming  through  a background 
of  field  particles  (labeled  31; 


slowing  down 


- v v ; 

s -a 


transverse  diffusion 


d / - n2  _ a/3  2 

! v - v ) = v v 

v a a x j.  a 


parallel  diffusion 


d / - \2  a/3  2 . 

— (v  - v ' = v v^  ; 

a 01  M H a 


energy  loss 


- v v. 


where  the  averages  are  performed  over  an  ensemble  of  test  particles  and 
a Maxwellian  field  particle  distribution. 


The  exact  formulas  may  be  written 


a/e 


= (1  + mc/n^)’!r(x  ) 


a/3x  a/e  -i 


v sec 

o 


a/p  r ./  a/e  w.  i /o  a/s.  , , ,,  ai/e>-i  a/e  -i 

= 2 u.tvx  )(1  - l/2x  ) + rl*  )J  v sec  ; 


a/e  _ r ,/  a/ew  a/e -i  a/e  -i  . 

= . inx  )/x  _ v sec  ; 

II  o 


a/e  _ 0 r / , s,,  a/e n 

Vg  = 2 L(  ma/tiu)i!i(x  ) 


,/  a/eN-i  a/e  -i 

I (x  ' " ®“<~ 


v sec  , 


where 


a/S  = Ve2e2>  n /m2^  sec"1 

0 a s ae  s a a sec 


x:i/3  = - nuv^/kT  ; 

~ p a e 


x 

ijf(x)  = ^ / dt>J~t  e C 

S11  J 


'!r7(x)  = dir/dx  , 


and  \ = is  the  Coulomb  logarithm  (see  below).  The  limicrng  forms 

QCu 

of  v . v and  v are  given  in  the  following  table.  The  two  expressions 
s i II 

Ci/3 

shown  for  each  rate  hold  for  very  slow  (x  « 1)  and  very  fast 

cy.  / s 

(x"  » 1)  test  particles,  respectively.  Test  particle  energy  e and 

field  particle  temperature  T are  both  in  eV;  a = m./m  , where  m is  the 

i p p 

proton  mass;  for  electron- electron  and  ion-ion  encounters,  field  particle 
quantities  are  distinguished  by  a prime. 


Electron- electron 


s/e’ 


/n  ,X  , 
e ee 


X 10 


-3/2 


7.7  X 10' 


-3/s 


ve/e'/n  ,X  , « 5.8  x 10" 3 T_1/2  e'1  - 7.7  X 10"3  e'3/2 

i e ee 

ve/e'/n  ,X  , «2.9  X 10"3  T-1^2  e"1  -3-9  X 10"3  Te'S/2 
11  e ee 


Electron- ion 


ve/l/n.Z2\  . 5s  0.2;  U3/2t'3/2  -3.9  X 10"6  e'"/2 
s 1 ei 


e/i . _s,  o c , .-/ 2~  x/2  -1 

v / n . Z X . « 2 . 5 x 10  u T S - 

1 1 te 


’.7  x 10"3  s'3/s 


£»/l  — *4^/2— l/->— —Q—  1 

v /n.Z‘X.  ss  1.2  x 10  u T /2e  - 2.1  x 10  u 'Te 

1 te 


Ion-electron 

vl/e/n  X.  « l.o  X 10"9  U_1T"3/2  - 1.7  X 10"4  u1/c  e'3/ 

s e le 


i/e 

v 


/n  \ . 

e le 


3.2  x 10" 


1.8  X lo'7  if 


1/ e 

v 

I! 


/n  1 . 

e le 


1.6  X 10" 


l.T  x 10" 


Ion- ion 


i/i ' 


vs  /ni'Z  Z'~\i'  *5.8x10- 


8 U 


r/: 


1 + 


T) 


9.0  X 10"8  u'l/s  (i 


-) 

u'  / 


i/ 1 f — o 

v /n  ,Z2Z'  x. . , 

i l li 


1.9  x IQ'7 


i -l/p  -3/- 

1.5  x 10  u 2e  /<i; 


v'71  /ni'z2z’\i'  *°'-8  X 10- 3 u^W^V1 


9.c  x 1C-3  Li1/2u'_1Ts'5/" 


In  the  same  limits,  the  energy  transfer  rate  follows  from  the  identi 


tv 


v£  - 2 vs  - vA  - v , 


except  for  the  case  of  fast  electrons  or  fast  ions  scattered  by  ions, 
where  the  leading  terms  cancel.  Here  the  appropriate  forms  are 

e/i 


v - L.P 

e 


.2  X 10- 


- S.Q  X 104  { u/T^  1 'Z  ~~1 


e exp(-  I*  us/T 


, 1 -1 
sec 


and 


i/i 


— 1.6  X 10 


n.Z2Z'  X.  . , 

i ii 


*3/2.  1 / : 


- 1.1  (u'/T)  1 exp( - u'e/TlJ  sec'*  . 

Cc/3 

In  general,  the  energy  transfer  rate  vg  is  positive  for  e > and 
negative  for  s < e^*,  where  xtt  ■ (m^/m.^  e */Tg  is  the  solution  of 
m^/rn^,  = 'i( x*') 'i 1 ( x^  . 

The  ratio  e^*/Tg  is  given  for  a number  of  specific  CX,  5 in  the  following 
table: 


Fig.  2 - Test  particle  showing-down  rates  fsa^,  scaled  by  the  ion-electron 
rate  esl'p(0)  in  the  limit  of  zero  relative  velocity  v,.  as  a function 
of  vt  divided  by  the  electron  thermal  velocity. 


The  rma 1 equilibration 


if  A < 0,  tan  is  replaced  by  tanh  " (-  A)1^2/(-  A)1  ^ 


For  T « T = T, 

1 ll 


v®  = 3.2  x 1CT7  n AT2/2 


= 1.9  x icrs  n>.z s/^r,z 


-i 


Coulomb  logarithm 

For  test  particles  of  mass  m , charge  e = Z e,  scattering  off 

Of  Cc  Of 

field  particles  of  mass  m^,  charge  e^  = Z^e,  the  Coulomb  logarithm  is 

defined  as  \ = in  A = £n(r  /r  . ) . Here  r . is  the  larger  of 

max  mm  mm 

e,.e3/(m^3u2’)  and  ^*/(2  m,^u),  averaged  over  both  particle  velocity 

distributions,  where  nm,  = m.nu/(m.  + m, ' and  u = v„  - v, ; 

Cep  up  u P ~ — 0t  — P 

r 3 (UrrEn  e~  /kT  ) where  the  summation  extends  over  all  species 

max  V Y V 

_ „ . l /- 

v for  which  u2  < v~  , with  v = (kX,/m  ) ^ . If  this  inequality  cannot 

Ty  TV  ' V 

be  satisfied  or  if  either  hot  “1  <r  or  u uu  J'  < r , the  theory 

cCc  max  cP  max 

breaks  down.  Typically  X « 10  - 20.  Corrections  to  the  transport 
coefficients  are  0(X-1),  hence  the  theory  is  good  only  to  ~ lof-  and  fails 


when  X ~ 1. 


The  following  cases  are  of  particular  interest 
(a)  Thermal  electron-electron  collisions 

X = 23  - in  (n  1/2  T 


T £ 10  eV 

e 


- 2k  - In  (n  1 /2  T -1' 

e e 


T * 1C  eV 
e 


4 


'b'  Electron- ion  collisions 


\ . = =2 
ei 


5 - in  ( n 1/2  ZT  ~3/2) 


> 10  eV  > T >T.m/m 

e i e 


2h  - in  (n  1/2  T ~ 1 ) , T > 10  eV  > T.m  /m.  ; 

<=  e e l e l 


' n;  1 Tr~^  zV1),  T.  > T m./m 


l l 


i e i e 


(c)  Mixed  ion- ion  collisions 


V'  = Vi = 25  ■ ** 


^ 1 2. 

"zz'(u  + u')  / nj_zS  n./z'2  \ | 

uT. ' + u'T.  \ T.  T.  ' / J 

*-  i i ' t l 


(d'  Counterstreaming  ions  (relative  velocity  vD  = 3 c)  in  the 
presence  of  warm  electrons,  klWn^  > v2  > kl^/m.,  kT^/m., 

V 


Xii'  = Vi  = 35  - in 


1 IC' 


zz'fu  + u'l 


uu'32 


(?) 


l 


Examples 

1.  Find  the  slowing- down  and  energy- trans fer  rates  for  a 100  keV 
deuteron  beam  in  a plasma  consisting  of  1 keV  electrons  and  100  eV 
tritons,  each  having  number  density  1014  cm 


Answer:  The  Coulomb  logarithms  are 


= 24  - In 
De 


[do14)1' 


from  formula  (b ' , and 


= 35  - lr 


(5/6^ 

(2.2  X 10£/3  x 1010)2 


from  formula  (dl.  The  beam  velocity  v^  = 2.2  x 10“  cm/sec  satisfies 

v_  « v_  « v , where  v_  and  v are  the  respective  thermal  velocities. 
L D e T e 

D/e 

It  is  therefore  appropriate  to  use  the  first  form  for  and  the 

, - D/T  „ 

second  for  v . Hence 


vD/e  = (1.6  X 10"9 )(2''1(l03'l"”/2(  1C14H15'> 


-1 

5 7 sec  , 


= (9  x 10-^(2f1/2U  + 2/3')(  10sr,3/2(  1014') ( lS'i 


3.9  sec 


The  net  rate  of  transfer  of  energy  to  the  electrons  is 


from  (b'1  and  (c),  respectively.  The  thermal  equilibration  rates  are 


thus 


vH  /e  = 3.2  X 10'9  n\  , T"3/s 

H+e  6 


2 -1 


4.5  X 10  sec  , 


and 


-H+/Li+++ 

v 


1.8  X 10' 


-19 


(7^1/2(32)(lO"2n)X 

H Li 


(1.7  X 10) 


;5  + ( 7'1  ( 100 'I  j 3 /S 


Hence 


= . 12  sec  “ . 


dT 


.+++ 


db 


H+  - H /Li 
- v 

e 


r i 

. H /e  r 

[25  - 100  J 

+ vg  500  - 100 

- 9 + 18  = 9 eV/sec  > 0. 


dT 


Evidently  as  T -*  <=  the  electron  contribution  to  — ; — vanishes  ~ T 

dt 

dv 

and  for  sufficiently  large  T < 0.  Thus  with  T =5  keV, 

6 dt 


-H  /e  , .3 

v = 1.4  x 10  sec 
e 


and 


dT 


dt 


( -12> 


25  - 100 J + (1>  x 10'°J)  5000  - loo] 


2 eV/sec. 


Hence  heating  the  electrons  cools  the  H ions! 


Fokker-Planck  equation 


+v.7^  + 

Dt  at 


f *7  r 

v 


-(£) 


coll 


where  F is  an  external  force  field.  The  general  form  of  the  collision 
integral  is( 


<*) 


coll 


= - y 7 .ja/3  with 

t ? ' 


2 2 

a/3  = :<S 

J - ^TrXC)p 

ma 


da  v r ( u‘ 


I - u u )u" 


J 1.  f^v)  V^f3  (v#)  - — f0(v#)7v 
”3  " Y ma  ” ‘ 


(landau  form)  where  u = v - v and  I is  the  unit  dyad,  or  alternatively 


a/3  , , 

J = 


eae3 


| f^v)  7 h(  v"1  - - 7 * [*  f ~^(  v)7  7 G(v)l, 
I • V - u V V 


where  the  Rosenbluth  potentials  are 


G(v)  =J'  f^(v/)ud~v/ 


H(v) 


(' 


u” 1 d' J v 


/ 


If  species  Ct  is  a weak  beam  (number  and  energy  density  small  compared 
with  background)  streaming  through  a Maxwellian  plasma,  then 


Ta/e  a/e  jx  l a/6  2_  jx 

j = - v vf  - — v v v r 


1 r a/e  a/e 


+ ± (v  - V 

2 1 II 


v . 7 f1"* 


Transport  coefficients ' - 

Transport  equations  for  a multispecies  plasma: 


a 

D ^ 


Dt 


+ na  -'-a 


1 


where  D '/Dt  = 9/ot  + v • 7 and  p.  -0 12 ) n^kT„  . R.v  = Z R,o  and  Q.  = J Q 

*V-v  s- v - - ^ v.V  “l  .O 

where  and  Q„g  are  respectively  the  momentum  and  energy  gained  by  the 

ath  species  through  collisions  with  the  is  the  stress  tensor,  and 

q,  is  the  heat  flow. 

Tiie  transport  coefficients  in  a simple  two-component  (.electrons  and 

singly  charged  ions'*  plasma  are  tabulated  below.  Here  and  i refer  to  the 

direction  of  the  magnetic  field  B = bB ; u = v^  - v^  is  the  relative 

streaming  velocity;  ng  = n^  = n;  j = - neu  is  the  current;  fjuce  = l.?t>  x 

10  B and  ou  . = ( m /m.'iau  are  the  electron  and  ion  gvrofrequencies , 
ci  e l ce 

respectively;  and  the  basic  collisional  times  are  taken  to  be 


N me(kTe^/2 

— . : = 3 .Mf  x 10 

b \ 2rr  n>  e^ 


T 

I ,13  ® 


3/s 


and 


3 \!m.(kT.  V 

< iv  l 

b J-  n>.e‘ 


= 2.09  x 107 


T. 


n\ 


•3/p 


i/s 


n). 


In  the  limit  of  large  fields  ((UcjTj  >:>  1'  the  transport  processes  may  be 
summarized  as  follows: 


momentum  transfer  R . = - R.  = R = R + ; 

-ei  -ie  - -u  -r 


frictional  force  R = ne(  j /cr  + j /r  N 

-u  -11  - |1 


ne 


conduct ivi t ies 


5.00  = 


r 


thermal  force 


ion  heating 


R = - 0.71nk7  T - - — — b x 7 T 
-T  -lie  p -13  T - -±  e 

- ce  e 

m nk 

Q • 3 3 — (T  - T ) ; 

i e i 

m.  t 
l e 


electron  heating  Q = - Q.  - R*u  ; 

e l - - 


ion  heat  flux 


ion  thermal 
conductivities 


q.  = - h17  kT.-  hL7  kT.  + hSj  x 7 kT . ; 

-i  11-11  i l-l  l a-  'll 


nkT.  t. 
h1  = 3 .9  — — 
m. 


nkT. 

l l 

H = 2 

i 


m.o)  .t. 

1 Cl  1 


5 nkT. 


A 2 m.uu  . 

1 ci 


e e e 

electron  heat  flux  q = q + q„,  ; 

- u - i 


3 nkT 


frictional  heat 
flux 


q = 0.71  nkT  u + - 
Hi  e-  i 


b x u ; 


X T 

ce  e 


q®  = - h6?  kT  - k*VW  - HCb  x 7 kT  ; 
-T  it -II  e i-i  e a-  -l  e 


thermal  gradient 
heat  flux 


electron  thermal 
conductivities 


stress  tensor 
(both  species) 


(here  the  z axis 


nkr  t 

*e  = 3.2  -£-S 

m 


nkT 

e _ ^ e 

K -4.7  


m 0)  t 

e ce  e 


c nkT 

e _ 5 e_ 

“ m a) 
'2  e ce 


P = (w  + w ) — - (w  - w ) - r w 

xx  p xx  yy  ? xx  yy  3 xy 


p = 2 (w  + W ) + (W  - W > 

yy  _ xx  yy  xx  yy 


n w ; 

2 XV 


P =P  = - Tl  W + -2.  ( w - W ' ; 


xv  yx 


l xy  2 xx  yy 


P = P 
xz  zx 


■ W - Tl  W 

2 XZ  4 VZ 


p = p * - 7 W + 7 W 
yz  zy  2 yz  4 xz 


P * - T|  W 
zz  o zz 

is  defined  parallel  to  BN  ; 


For  both  species  the  rate-of- strain  tensor  is  defined  as 


cv . Sv,  2 

W = — 1 + — - - - 5.,  7-v 

Jk  5 


When  B = 0 the  following  simplifications  occur: 


R » nej/c  ; 
-u  - 11 


R = - 0.71  n7kT  ; 
-X  e 


q.  = - ;i  VkT.  ; 
-i  ii  - i 


qu  = 0.71  r*XeU  J 


qfi  = - >te  7 kT  ; 

~T  II  - e 


P..  = - 7 W..  , 

j k o jk 


when  id  t » 1 » uu  .t.,  the  high-field  expressions  are  obeved  b'  the 
ce  e ct  l 

electrons  and  the  zero- field  expressions  by  the  ions. 

Collisional  transport  theory  is  applicable  when  (l)  macroscopic 
time  rates  of  change  satisfy  ^ « — , where  - is  the  longest  colli- 
sional time  scale,  and  fin  the  absence  of  a magnetic  field'  (C'l 
macroscopic  length  scales  L satisfy  L » £,  where  l = v-  is  the  mean 


l 

L 


ft 


free  path.  In  a strong  field,  t»  1,  condition  (2*1  is  replaced  by 

L » l and  L »\£r„  (L  » r in  a uniform  field'',  where  L is  a 

macroscopic  scale  parallel  to  the  field  B and  L is  the  smaller  of 

(V^B/B^“~  and  the  transverse  plasma  dimension.  In  addition,  the  standard 

transport  coefficients  are  valid  only  when  (3)  the  Coulomb  logarithm 

satisfies  X » 1;  (^'  the  electron  gyroradius  satisfies  r^  » X , or 

B2  « Srrn  m c~  ; (5^  relative  drifts  u = v - v„  between  two  species  are 
e e -a  -3  r 

small  compared  with  the  thermal  velocities,  u «(kT  /m  r/2,(kT,/nOl/2, 
and  (o'  anomalous  transport  processes  owing  to  microinstabilities  are 
negligible . 


Weakly  Ionized  Plasmas 

Collision  frequency  for  scattering  of  charged  particles  by  neutrals 


n c (kX  /m 

o OCo  3 a 


l/„ 


where  n is  the  neutral  density  and  c is  the  cross-section,  typically 
o jCo 


~ 5 x 1C"1-  cm"  and  weakly  dependent  on  temperature. 


i 


i 

■j 


i 


T D +T  D. 
i e e l 


cm  /sec 


where  - e^/m^y^  is  the  mobility.  The  conductivity  c.  satisfies 

Ca  ~ naeolJ'a  ' 


In  the  presence  of  a magnetic  field  B,  U and  - become  tensors 


ra  0!  _ 

J = c *E 


O^E  + caE  + a" E x b , 

ll-ll  A* 


where  b = B/B  and 


a 

II 


naea/nava 


a 


a 2//  2 
c v A vf 
II  ^ 0! 


,.2  'I 


V it  / ( v2 

a ca  v a 


+ -i)2  ) 

ca 


Here  c.  an<^  p are  t^le  Hall  and  Pedersen  conductivities,  respectivelv 


IONOSPHERIC  PARAMETERS 


These  are  typical  nighttime  values.  Where  two  numbers  are  entered,  the 
first  refers  to  the  lower  and  the  second  to  the  upper  portion  of  the  layer. 


Quancity 

E Region 

F Region 

Number  density  (m"3) 

1.5  X 1010  - 3.0  X 1010 

5 x 1010  - 2 x 1011 

Height-integrated  number 

9 x lo14 

4.5  x 1015 

density  (ra"2) 

Ion- neutral  collision 

2 X icr  - 10s 

0.5  - .05 

frequency  (sec”1) 

Ion  gyro-/collision 

.09  - 2.0 

4.5  x 102  - 5.0  x 103 

frequency  ratio  v. . 

Ion  Pedersen  conductance 

0.09  - o.5 

2.2  x 1CT3  - 2 x 1CT4 

*./(l  + *?  ) 

1 1 

Ion  Hall  conductance 

£ X icr4  - o.a 

1.0 

1 + 

l l 

Electron- neutral  collision 

1.5  x 1C4  - 9.0  x 102 

SO  - 10 

frequency 

| Electron  gyro-/collision 

4.1  X 102  - 6.9  X lo3 

7.8  x 104  - 0.2  x 105 

frequency  ratio  a 

e 

1 Electron  Pedersen 

2.7  X lCT3  - 1.5  X 1C*4 

1J-5  - 1.5  X 10- 

conductance  •'  /(  1 + ) 

e e 

t Electron  Hall 

1.0 

1.0 

conductance  aA/(  1 + ^ 

e e 

Mean  molecular  weight 

28  - 26 

22  - 16 

ton  gyrof requencv  sec"* 

1-0  - 190 

2J0  - l00 

| Neutral  diffusion 

T0  - X let 

io'- 

1 coefficient  nr/sec'' 

J 



The  terrestrial  magnetic  field  in  the  lower  ionosphere  at  equatorial  lati- 
tudes is  approximately  • 1'"  tesla.  The  earth's  radius  is 

q « - v IT/^  km  : > I ' kn . 

E 


4 


1 


fci 


4 


CTR 


L20J 


Natural  abundance  of  deuterium  n^/n^  = 1-5  x 10  4 
Mass  ratios 

me/mD  * 2 -72  x 10"4  = 1/3o70 

(m^/nr  = 1.65  x 10  2 = l/oO.o 

m /ittj.  = 1.82  X 10‘4  = 1/5496 

6 i-  P 

(me/mT)  = 1.35  x 10"2  = 1/74.1 


Fusion  reactions  (branching  ratios  are  correct  for  energies  near  the 

. [2  2] 

cross-section  peaks;  a negative  yield  means  the  reaction  is  endothermic). 


( la) 

(lb) 

(2) 

(3) 

(4) 

( 5a) 

(5b) 

(5c) 

(6) 

(7a^ 

(7b  'l 

(8) 

(9) 


D + D 


50$ 

50! 


+ T 
+ He3 
+ T 


D 
D 
T 

He3  + T 


+ Li'3 
+ Li7 


51$ 


43$ 


6$ 


20" 


Li': 


T(1.01  MeV)  + p(3 .02  MeV"* 

He3  ( . S2  MeV 'i  + n(2.45  MeV) 
He4(3 • 5 MeV)  + n(l4.1  MeV) 
He4(3.o  MeV)  + p(l4.7  MeV) 
He4  + 2 n + 11.3  MeV 
He4  + p + n + 12 . 1 MeV 

He4(4.S  MeV)  + D(9.5  MeV) 

He5(  2.4  MeV 'I  + p(ll.9  MeV) 

He4( 1.7  MeV)  + He3(2.3  MeVl 
2 He4  + 17.3  MeV 


Be'  + n - 1.6  MeV 


~ 50$ 


+ B 


11 


2 He4  + 22.4  MeV 

3 lie4  + -.7  MeV 


uming  the  target  ion  at  rest,  can  be  fitted  by 


■4 


' 


For  low  energies  (T  ^ 25  keV)  the  data  may  be  represented  by 

2 / 1 j 

(tv)pD  = 2.53  x 10  T 3 exp(  - 18.76  T 3)  cm3/sec 

(cv)q^,  .=  3-68  x 1CT12  T exp(  - 19.9^  T ) cm3/sec 

The  power  density  released  in  the  form  of  charged  particles  is 

= 3-3  x 10  1 n^Cov)  watt/cra3  (including  subsequent  D-T  reaction) 

PQT  = 5.0  x 1CT13  n^Cov)  watc/cm3 

?DHe3=  2‘9  X 1Cr13  W.e^W  watc/cm' 

The  curie  (abbreviated  Ci)  is  a measure  of  radioactivity:  1 curie  = 

3.7  x 10^°  counts/sec.  Absorbed  radiation  dose  is  measured  in  rads: 

1 rad  = 102  erg/g. 


— . 


RELATIVISTIC  ELECTRON  BEAMS 


Here  y = (l  - B2)"^  is  Che  relativistic  scaling  factor;  in  analytic 
formulas  units  are  mks  or  cgs,  as  indicated;  in  numerical  formulas,  I 
is  in  amps,  B in  gauss,  electron  density  N in  cm-1,  and  temperature, 
voltage  and  energy  in  MeV  ; B = v /c  . 

Z 2 


Relativistic  electron  gyroradius 


Relativistic  electron  energy 
Alfven- Lawson  limit 


r = f - 1)  E ( cgs) 

6 eB 

1 

= 1.70  x lCPtv2  - 1)*  B cm 
W = mc^y  (cgs)  = 0.511  Y MeV 

3 

_ inc  o / \ 

IA  * (cgs) 

e 

= — —By  (mks) 

u e 2 
o 


= 1.70  x 104  BzY  amp 


The  ratio  of  net  current  to  I is 

I/IA  = v/y  , 

where  v = Nr  , with  r = e‘/mc'  = 2.82  X 10  13  cm.  Beam  electron  number 
e e 

density  is 

n = 2.08  X 10 ? J/B  cm*’, 


where  J is  current  density  in  amp/cm  . 


For  a uniform  beam  of  radius  a 


Child's  law  (space-charge  limited  current  density  between  parallel 
plates  with  voltage  drop  V in  MeV  and  separation  d in  cml 


Nonrelativist ic : 


J = 2.3^  x 1C P V /2  d''r 


Relativistic  (y^  2j; 


J = 2.71  X 103  (vlA 


>.ci*7l)2  d'C  ^2 


Volf-’ge  registered  by  Rogowski  coil  of  minor  cross-sectional  area  A,  n 
turns,  major  radius  a,  inductance  L,  external  resistance  R and  capacitance 
C all  in  mks''  : 


self- integrating 


R nAU  I 

V volts 

L 2na 


1 nAu  I 
o 

externally  integrated  V = — volts 

RC  2TTa 


X-ray  production  for  target  with  average  atomic  number  Z ( V ^ 5 MeV) 


7 = x-ra?  P°wer  = 7 x icr4  zv 

beam  power 


X-ray  dose  at  1 meter  generated  by  an  e-beam  depositing  total  charge  Q 

coulombs  while  V > .54  V in  material  with  charge  state  Z 

max 


D = 150  V'  ’ Q Z 

max 


1 


ELECTROSTATIC  STREAMING  INSTABILITIES 


[24] 


In  this  table,  subscripts  e,i,d,b  stand  for  "electron,"  "ion,"  "drift," 
and  "beam,"  respectively.  Thermal  velocities  are  denoted  by  a bar.  In 
addition,  the  following  are  used: 


I 


m 

M 

V 

T 

n , n. 
e 1 


C = T /M 
s e 


U)  , U>. 

e t 


r , r. 
e l 


8 

V. 


W 9 ... 

e 1 


*“H 

U 


electron  mass 
ion  mass 
velocity 
temperature 
number  density 

harmonic  number 

ion  sound  speed 

plasma  frequency 

Debye  length 

gyroradius 

plasma/magnetic  energy  density  ratio 
Alfven  speed 
gyro  frequency 

hybrid  gyro  frequency,  = .T  L 

H e 1 

relative  drift  velocity  of  two  ion  species 


I 


Parameters  of  Most  Unstable  Mode 


LASERS 


Parameters  of  Important  Laser  Systems 

Efficiencies  and  power  levels  are  approximately  state-of-the-art  ( 1976) 


Type 

Wavelength 

(microns) 

Efficiency 

Power  level 

s available(Wl 

Pulsed 

CTv 

CO 

2 

10.6 

.01  - .02  pulsed 

> 1 o11 

> 103 

CO 

5 

oA 

> 10? 

> 100 

Iodine 

1.515 

3 x 10-3 

> lO1- 

— 

Ud-glass,  YAG 

1.06 

10"3 

> 1012 

1 - 300 

Ruby 

.69^3 

< 10‘3 

l'O10 

1 

He  - Ne 

. o32  c 

10'4 

— 

1 - 50  X 10~3 

Argon  ion 

M - .60 

io"3 

5 X 104 

1 - 10 

N 

.3371 

lo"3  - .05 

»— * 

0 
01 

1 

►— * 

O 

2 

Xenon 

.175 

.00 

> 10a 

5.  Ponderomotive  force 


where  ± refers  to  the  velocity.  The  rate  of  change  of  polarization 
angle  S as  a function  of  displacement  s (Faraday  rotation)  is  given  by 


— = - k(n  - n ) - 2.55  NBf'2  cm'1  , 
ds  2 

where  N is  the  electron  number  density,  B is  the  field  strength,  and  f 
is  the  wave  frequency.  For  N = 1C^°  cm  , B = 10  kG  and  f 3 5 X 1.'  Hz 

—3 

(corresponding  to  a wavelength  \ = 1 u),  — - - 2 .0  x 10  rad/m. 

. ds 

2.  The  quiver  velocity  of  an  electron  in  an  e-m  field  of  angular 
frequency  jq  is 


eE  _ 

v = — = 25.6  \1  I \ cm/sec 
o ’ o 


ATOMIC  AND  RADIATION  PROCESSES 


The  subject  ot  atomic  physics  is  central  and  indispensible  for 
unstanding  the  behavior  of  real  plasmas.  It  is  ill-suited  for  dis- 
tillation into  a few  general  formulas,  however.  Those  compiled  here 
should  not  be  applied  blindly.  The  inexperienced  user  would  do  well  to 
first  consult  the  references  cited,  particularly  _31  _ and  _3 5 which 
are  more  accessible  to  plasma  physicists. 

Energies  and  temperatures  are  in  eV;  all  other  units  are  cgs 
except  where  noted.  Z is  the  charge  stat • , Z =0  refers  to  a neutral 
atonO;  the  subscript  e labels  electrons.  N refers  to  number  density, 
n to  principal  quantum  number.  Asterisk  superscripts  on  level  popu- 
lation densities  denote  local  thermodynamic  equilibrium  LTE ^ values. 
Thus  N^  is  the  LTE  number  density  of  atoms  or  ions'  in  level  n. 


Characteristic  atomic  collision  cross  section 

1 rra2  = 3. SO  x 10-17  cm2 

o 

Binding  energy  for  outer  electron  in  level  labelled  by  quantum  numbers 

n,  H 

z2eh 

(2)  (M)  = - ■ eV, 


where  E^  = 13.6  eV  is  the  hydrogen  ionization  energy  and 


l > 5,  is  the  quantum  defect. 


Excitation  and  Decay 

Cross  section  (Bethe  approximation)  for  electron  excitation  by 

[25,2.'] 

dipole  allowed  transition  m -♦  n ‘ 


(3) 


a = 2.56  x 10"13 
mn 


/ g(n,m 

r-i  ' 


AE 


cm- 


tun 


where  / is  the  oscillator  strength,  g(n,m)  is  the  Gaunt  factor, 
run 

€ is  the  incident  electron  energy,  and  AE  = E - E . 

run  n m 

Electron  excitation  rate  averaged  over  Maxwellian  velocity  distribution 


l2T,2-J 


(M 


mn 


1.6  x 1 O'5/  <g(n,ir.'>N 

N < c v > = — -exp 

6 11111  ^E  vT" 

run  e 


/ -E  \ 

(-  —jsec  x. 


where  < g(n,m)>  denotes  the  thermal  averaged  Gaunt  factor  whose  behavior 
•as  a function  of  energy  is  shown  in  Fig.  3 (generally  ~ 1 for  atoms, 

~ 0.2  for  ions) . 

Rate  for  electron  collisional  de-excitation 


(5) 


Y = (N*/N*)X  sec'1, 
ran  tn  n mn 


Here  /N,t'  = g /g  'exp'AE  /T  ) is  the  Boltzmann  relation  for  level 
m n m n run  e 

population  densities,  where  g^  is  the  statistical  weight  of  level  n. 
Rate  for  spontaneous  radiative  decay  n *♦  m (Einstein  A coefficient  ^ 


A * 4.*>  x 10  7 — I u-E  sec 
iin  " & run  run 

m 


1 


Fig.  3 — The  effective  Gaunt  factor  g for  neutral  atoms  and 
positive  ions.  For  neutrals  the  curve  is  deduced  from 
the  experimental  results  for  Is  —■  2p  in  H. 


Intensity  emitted  per  unit  volume  from  the  transition  n m in  an 
optically  thin  plasma 


I = 1 ,6  x 10~19A  H 4E  w/cr 
run  run  n nm 


Condition  for  steady  state  in  a corona  model 


NN  < c v > = N A , 
o e on  n no 


where  the  ground  state  is  labelled  by  a subscript  zero. 

Hence  for  a transition  n - ji  in  ions,  where  < g(n,m')  > « 0.2, 


(9)  Im  - 5-1  x IQ'25  fnra(80/gm  (AEnm/AEno)3NeNoTe-1  2 exp(- AE^/l^  W/ 


Ionization  and  Recombination 

In  a general  time- dependent  situation  the  number  density  of  the 
charge  state  Z satisfies 


= N [ - S'Z)N(Z)  - cv ■'  Z )N  ( Z 4 + S^Z  - l)N(Z  - 1' 


+ ar(Z  + 1 )N  ( Z + 1)3. 


Here  S^Z)  is  the  ionization  rate.  The  recombination  rate  avZ)  has  the 
form  or(Z)  = a (Z)  + Ng  0/3(2),  where  and  23  are  the  radiative  and 
three-body  recombination  rates,  respectively. 


Classical  ionization  cross-section  for  any  atomic  shell  k 


. 1n-,,  \ 

Oj  = o x 10  jp  gk  X cm 
k 


■1 


‘ ! ■'  m 


- ■ — — 7 — ^ 


1 " ' 


— - 





Here  is  the  number  of  shell  electrons,  is  the  binding  energy 
of  the  ejected  electron,  x = e/U^  where  e is  the  incident  electron 

energy,  and  g is  a universal  function  with  a maximum  value  ss  0.2  at 

1 

X « 4. 


Ionization  from  ion  ground  state,  averaged  over  Maxwellian  electron 
[*2  8 1 

distribution  ' 


10'S(T  /E^)1/2  z 

(12)  S (Z  ) = -S- exp  (-  E /T  ) cm3/sec, 

(E b3/2(6  + T /E I)  6 

v CO  ' 0 CO  ' 


where  E^  is  the  ionization  energy. 

Electron-ion  radiative  recombination  rate  (e  + N(Z)  N(Z-l)  + hv ) -30- 


(13) 


a ( Z)  = 5.2  x 10~14  Z JeZ/T  \o.k-3  + 

r 00  6 I 

4 4n  (EZ/T  ) + 0.k69  (EZ/T  )_1/3j  cm3/sec. 
2 00  e 08  e i 


|2  0 I 

For  10  keV  < Tg/Z2  < 150  keV,  this  becomes  approximately1 

«■  72 

(14)  ot  (Z)  = 2.7  x lO*”13  cm3/ sec . 

r ^ 

13  2 | 

Collisional  (three-body'  recombination  rate  for  singly  ionized  plasma1 


* 15)  ^3  = 8.75  x lO-27  T ~9^  cm  /sec 

e 

limited  to  low  temperatures,  kT  ■'  e"  ^ • Photo ionizat ion  cross 

0 CO 


. 


--  — — 


section  for  ions  in  level  n,  SL  (short  wavelength  limit') 


( lo)  c ^(n,  V)  = 1.64  x 10  Z“/(n3Kr  +~Z)  t 

where  K is  wavenumber  in  Rydberg  units  (l  Rydberg  = 1.0974  X 105  cm  X). 
Charge  Exchange 

Two  kinds  of  charge  exchange  take  place  in  a plasma,  depending  on 

AE  in 

(17)  A + B+  - A+  + B + AE  . 


When  AE  = 0 (the  internal  energy  difference  vanishes^  the  process  is 
called  resonant.  When  A = B the  resonance  is  symmetric;  otherwise  it 
is  asymmetric.  When  AE  f4  0 the  process  (17)  is  just  a charge  transfer 
reaction.  In  this  case  the  adiabatic  principle  of  Massey  suggests  that 
the  peak  of  the  cross  section  c arises  when 


V 


Here  £ is  the  interaction  range  and  V is  a critical  relative  velocity. 
For  larger  values  of  V, 


(18) 


k exp 


/ el 'e! 

\ 4 hV 


with  k being  a constant  for  each  reaction.  However,  for  symmetric 
charge  exchange  the  cross-section  as  given  by  Fursov  can  be  expressed 


32 


for 


Values  of  a and  b are  given  below  for  some  selected  reactions 

the  energy  range  E . < E < E 

mm  max 


Process 

a 

( 10  ' cm) 

b 

( 10  ~ cm) 

E . 
mm 

(eV) 

E 

max 

(eV) 

H+  — H 

7-6 

1.06 

20 

10,000 

He+  — He 

1.1 

20 

6,000 

Ne+  — Ne 

6 . 1 

1.7 

20 

6,000 

0 

1 

+a: 

• k.2 

O.36 

25 

£,000 

Equilibrium  Models 

Saha  equilibrium 


[33] 


(20: 


N N (Z) 
e 1 

N*(Z-1) 

n 


6.0  x 1021 


gZT  3/2 

si  e 

2-1 


exp 


cm 


7 

where  is  the  statistical  weight  for  level  n of  charge  state  Z and 

7 

E»(n>^)  is  ionization  energy  of  the  neutral  atom  initially  in 
level  (n,-0,  given  by  Eq.  (2). 


In  a steady  state  at  high  electron  density 


(21 


N N*(Z)  S(Z-l) 

e _ 

N*(Z-1)  = *3 


a function  only  of  T. 


gk 


Conditions  for  LTE 


(a)  Collisional  and  radiative  excitation  rates  for  a level  n must 

[33] 

satisfy 


(22  ) Y/A  nm  > 10  . 

tun  run 


(b)  Electron  density  must  satisfy 


^ 7 x 1013Z 


7n-17/2(T/EZ)1/2  cm-3 


Steady  state  condition  in  corona  model 


N(Z-l)  _ r 
'n(Z)  S(Z-l) 


Corona  model  is  applicable  if 


1012  t'1  < N < 1016T  7/2  cm-3  , 
I e e 


where  t'1  is  the  inverse  ionization  time. 


Radiation 

Average  radiative  decay  rate  of  state  with  principal  quantum  number  n is 


A - Z A = 1 .6  x 1010  Z'Si  9//s  sec”1 
n m < n nm 


Natural  line  width  (AE  in  eV ) 


AE  At  * h = 4.1^  x 10-1  5 eV-sec. 


where  At  is  the  lifetime  of  the  line. 


Doppler  width 


i-t/i  _ _ _ _ - 

7*  = 7 «7  x 10 


# - 


where  p,  is  the  emitting  atom  or  ion  mass  in  units  of  the  proton  mass. 


Optical  depth  for  a Doppler  broadened  line 


T = I.76  x 10"13  X = 5.4  x 10”9  \(j±y/2NX, 


where  \ is  wavelength  and  l the  physical  depth  of  the  plasma;  M,  N and 
T are  atomic  mass,  number  density  and  temperature  of  the  absorber;  p, 
is  M divided  by  the  proton  mass.  Optically  thin  means  t < l. 


Note  that  the  units  employed  here  are  those  summarized  on  p.S£. 


Resonance  absorption  cross  section  at  center  of  line 


(50)  = 5-6  x 10  ^ cm-  . 

c 

Wien  displacement  law:  wavelength  of  maximum  black  body  emission  is 

given  by 


> = 2.50  x lO-5  T_1  cm. 

max 


Radiation  from  surface  of  black  body  of  temperature  T 


W * 1,05,  X'  l55  T4  watts /cm2 


Bremsstrahlung  from  hydrogen- like  plasma 


P = 1.69  x 10"32  N T 1/2  S [Z2N(Z)1  watts/cm3, 


where  the  sum  is  over  all  ionization  states  Z. 


Bremsstrahlung  optical  depth' 


5.0  x 1 0~36  N NjZ2  g^t"7/2. 


where  g **  1.2  is  an  average  Gaunt  factor  and  l is  the  physical  path 


length. 
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I 

! 


Recombination  (free-bound)  radiation 


P = 1.69  x 10-32N  T ^2  E Z2fl 

r e e 


watts/cm3 , 


Cyclotron  radiation  in  magnetic  field  BL 


P * 6.21  x 10~21B2N  T watt/cm3 
c e e 


For  NgkTe  = N^kl^  = B2/16tt  (B  = 1,  isothermal  plasma). 


(37)  P = 5.00  x 10-31N  2T2  watt/cm3. 

' c e e 


Cyclotron  radiation  e-folding  time  for  a single  electron 


_ __  9.0  x 103b'2 
C “ 2.5  + V S6C 


..  r*V4 


where  y is  the  energy  dtv^|«Lby  me2. 


Number  of  cyclotron  harmonics  traj^^^tn  a medium  of  finite  depth 


n*  = (57  S Bi)l/a  , 


where  3 = NkT/SniJ2. 


Line  radiation  is  given  by  summing  Eq.  9 over  all  species  in  the  plasma. 
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